Femtosecond quasiparticle relaxation dynamics and probe polarization anisotropy in 

YSr^Ba^C^Og (a; = 0,0.4). 
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Femtosecond pump probe experiments are reported on quasiparticle relaxation and recombination 
in YSr2:Ba2-iCu408 as a function of temperature and polarization. The data show a 2-component 
relaxation similar to YBa2Cu307_i, one component being assocated with the superconducting tran- 
sition, and the other with the pseudogap below T* . The relaxation time r p associated with the 
pseudogap is found to be T-independent, while the relaxation time r 9 of the component obseved 
only below T c exhibits a clear divergence near T c . A strong polarisation anisotropy of the picosecond 
transient is observed below T c which is attributed to the anisotropy of the probe transition matrix 
elements. 



Femtosecond time-domain spectroscopy is capable of 
giving important information on the quasi-particle ex- 
citations and a low-energy structure of correlated elec- 
tron systems, and particularly high-T c superconduct- 
ing cuprates (HTSC). Measurements on YBa2Cu307_5 
(Y123) at different doping levels have shown two- 
component relaxation dynamics on the picosecond 
timescale, that was attributed to the simultaneous ex- 
istence of two gaps in optimally doped and overdoped 
regions Q. Moreover, similar two-component relaxation 
dynamics has been thus far observed in several HTSC 
H !> @> !> |> 0> §> suggesting this being a general fea- 
ture in these materials. The magnitude and sign of the 
two components at temperatures below T c were found to 
depend on the material, probe wavelength Q, and also 
the effects of probe polarization dependence were inves- 
tigated on untwinned Y123 single-crystals where a 
response parallel and perpendicular to the Cu-0 chains 
was separately probed and found to be different. In 
some cases (T12201[|, Bi2212 j|) the signs of the dif- 
ferent components observed in the relaxation were oppo- 
site. However, the anisotropy of the photoinduced signal 
with respect to the probe pulse polarization has not been 
discussed in detail thus far. 

The purpose of this paper is two fold. First, we show 
that the 2-component relaxation behavior observed in 
Y123 is similar also in YBa2Cu40s (Y124), and second, 
utilizing the fact that Y124 has a fixed oxygen content 
and a well-defined untwinned orthorhombic structure, to 
probe the polarization dependence of the photoinduced 
signal. 

We report on a systematic investigation of Y124 and 
YBai.6Sr .4Cu4Os (Y124:Sr) using the usual femtosec- 
ond time-resolved pump-probe technique. As discussed 
in detail previously Jll|], a short (~80fs) pump laser pulse 
excites the carriers in the sample. Photoexcited electrons 
and holes with energies on the order of photon energy 
quickly thermalize via electron-electron and electron- 
phonon thermalization, reaching states just above the 
gap in a time short compared to the pulse duration. The 
gap in the density of states presents a relaxation bot- 



tleneck, and the relaxation of photoexcited carrier den- 
sity near Ep is measured through measurement of small 
changes in the optical reflectivity AR/R or transmittance 
AT/T of the sample as a function of a time-delay be- 
tween the pump and probe pulses. In these experiments, 
a Thsapphire mode-locked laser, which operated at a 78 
MHz repetition rate, was used as a source of both pump 
and probe light pulses. The wavelength of the pulses 
was centered at approximately A « 800nm (1.58eV) and 
the intensity ratio of pump and probe pulses was ap- 
proximately 100:1. The pump and probe beams were 
crossed on the sample's surface, where the angle of in- 
cidence of both beams was less than 10°. The diameter 
of both beams on the surface was ^100/im and the sur- 
face was parallel to a ab-crystal plane. The typical en- 
ergy of pump pulses was 0.2nJ (1.25x 10 9 eV), which pro- 
duces a weak perturbation of the electronic system with 
approximately 3xl0 10 thermalized photoexcited carriers 
per pulse. (The approximation is based on the assump- 
tion that each photon creates Hu>/A thermalized pho- 
toexcited carriers, where A sa 40meV is of the order of 
the superconducting gap.) The train of the pump pulses 
was modulated at 200kHz with an acousto-optic mod- 
ulator and the small optical changes were resolved out 
of noise with the aid of phase-sensitive detection using 
EG&G digital lock-in amplifier model 7265. The pump 
and probe beams were also cross-polarized to reduce scat- 
tering of pump beam into the detector (avalanche photo- 
diode). A detailed description of the experimental tech- 
nique and the theory of excitation and relaxation of the 
photoexcited carriers in superconductors with different 
gap structures can be found in refs. |ll], 

The two Y124 samples used for this investigation were 
prepared in Zurich by a nonstoichiometric flux-growth 
technique, using a BaO-CuO eutectic mixture as the flux. 
In the Sr-doped compound, part of the Ba was substi- 
tuted with Sr. The crystals grew at pressure of 900bar 
and temperatures, which ranged from 1000°C to 1120°C. 
So obtained crystals had a thickness of approximately 
100/im in c-axis direction and the dimensions of 0.4x0.2 
mm 2 . The crystal axis were determined by x-ray anal- 
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ysis. Details of the growth method and the techniques 
used for characterization of the samples are given in ref. 

s 

The photoinduced (PI) reflection AR/R as a function 
of time at different temperatures is shown in Fig. 1 for 
the direction of polarization of probe pulse along a-axis 
and 6-axis of Y124. The chains are parallel to the 6-axis 
and the PI response in this directions clearly shows a 
presence of at least two relaxation process with differ- 
ent sign of AR/R. At temperatures T > T c we see a 
signal with a positive AR/R and with a relatively fast 
relaxation time (jp ~0.2ps). As the temperature is low- 
ered below T c , a second component with a longer re- 
laxation time {tq «2ps) and negative sign starts to ap- 
pear. This kind of behavior of the two components is 
also present in PI response with the probe polarization 
along a axis, which can be clearly seen, if data is pre- 
sented on a logarithmic scale as in the insert of Fig 1. 
The relaxation of AR/R after 200fs can be modeled with 
a function of the form AR(t,T)/R = G(T) exp(-t/r G ) + 
P(T) exp(— t/rp), where the temperature dependant am- 
plitudes G(T) (G(T) = for T > T c ) and P(T) can have 
the same or opposite sign below T c , depending on the 
polarization direction of the probe pulse. Fits, which are 
presented by a continuous line in Fig. 1, are in a good 
agreement with the data. We note that above T c the 
data shows a slight departure from a simple exponential 
decay. Similar behavior above T c has also been reported 
in ref. |^| , where two components of opposite sign below 
T c and one component slightly non-exponential above T c 
were observed on Tl 2 Ba2Ca2Cu30io- Whether a stretch 
exponential decay is relevant or the presence of an ad- 
ditional component is the reason for the discrepancy, we 
cannot assert with a good degree of certainty from the 
present data. 

In addition to the two picosecond components, a weak 
oscillatory component with frequency of 3.0(±0.3)THz 
can also be seen above T c . This is shown in Fig. 1 
c) . Such an oscillatory component of the same frequency 
was already observed with pump-probe measurements on 
YI24 at room temperature Oj. 

We analyze the temperature dependences of relaxation 
time tg and of both amplitudes G(T), P(T) in the same 
way as was done previously for YI23 Qj. The temper- 
ature dependance of quasi-particlc recombination time 
in a superconductor with a temperature dependant gap 
A(T) below T c is described by equation [O 



tg = 
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where oj p h is a typical phonon frequency, N(0) is the 
density of states (DOS), r w is a characteristic phonon 
linewidth and A(0) is the value of the gap at T=0 K. 
In Fig. 2 we show divergent behavior near T c in the 
temperature dependance of rg, predicted by the Eq. 
(1) (solid curve). The data points were obtained from 
the fits of the time evolution of AR/R for Y124 along 



the a and b axes and for Y124:Sr along the b axis. In 
the fits we have used the following values of parame- 
ters in Eq. (H) u ph = 400 cm" 1 , T w =10 cm" 1 , N(0) 
=5 eV _1 cell spin -1 and a BCS functional form for the 
temperature dependant gap A(T). The fitting parameter 
was A(0). 

In Figure 3 we compare the measured temperature de- 
pendence of the signal amplitudes G(T), P(T) with those 
predicted by the theoretical expressions p"2| for BCS tem- 
perature dependant gap 



G(T) 
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and for temperature independent gap 



P(T) = 
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exp[-A P /fc B T] 
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In the fits we used v =18 for the number of modes in- 
teracting with quasi-particles, f2 c =0.1eV for a typical 
phonon cutoff frequency and A^(0) =5 eV -1 cell -1 spin -1 
for the DOS. The agreement between the data and the 
theory is seen to be very good. 

In Fig 3 c) we show a polarization dependence of both 
amplitudes G(T, 6), P(T, 6) for Y124 at T = 45 K. P(T) 
does not show any polarization dependance, which was 
also confirmed by polarization measurements above T c 
In contrast, the polarization dependance of G(T) shows 
a significant angular dependence, exhibiting a change of 
sign as the polarization direction is changed from parallel 
to perpendicular to the Cu-0 chains. 

This angular dependance can be understood in terms 
of the anisotropy of the probe transition matrix elements. 
We first note that in spite of the fact that Y124 is or- 
thorhombic, the main contribution to the probe signal 
comes from interband resonance transitions [9] which in- 
volves atomic wave functions in which the dipolar matrix 
elements have a 4-fold rotational symmetry around the 
c axis. Without specifying exactly which transitions are 
involved, we can write the absorption coefficient in terms 
of the Fermi golden rule: 



a cx 



deN(e)N(e + uj)\M(e,uj)\ 2 f(e)(l - f(e + u)) 



Here N(e) is the density of electronic states, /(e) is the 
distribution function for holes, hio is the energy of the 
probe photons and M(e, uj) is the dipole matrix element 
for the transition. For small perturbations, AR is pro- 
portional to the photoinduced absorption Aa, so: 

AR cx Aa cx / deN(e)\M(e,uj)\ 2 (f' (e) - f(e)) (4) 



where / (e) is the nonequilibrium distribution function of 
the charge carriers. The integral is taken in the vicinity of 
the Fermi energy Ep over the width of the resonance [9] 
and we assumed for simplicity that N (e + f2) is constant 
within the resonance width. 
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If M(e, u>) is constant over the whole range of energies 
— TiLo < e — Ep < hoj, then AR = Aa = because of the 
conservation of particles. In other words, the photoin- 
duced increase in absorption due to probe transitions 1 
(see Figure 4) originating from the photoexcited electron 
states at e — Ep > A just above the gap to unoccupied 
hole states at e — Ep ~ hio + A exactly cancel the de- 
crease in absorption due to transitions 2 originating from 
occupied electronic states at e — Ep < — A just below the 
gap to unoccupied states at e — Ep ~ tuo — A. The same 
is true for the hole transitions 3 and 4 (see Fig. 4). 

Turning to the situation in hand, we assume that 
M.(e 7 uj) is not constant and proceed to derive an ex- 
pression for the polarization anisotropy of the probe ab- 
sorption. The general expression for the square of the 
dipolar matrix element can be written as: 

|M(e)| 2 = M 2 x (e) sin 2 (9) + Af 2 (e) cos 2 (9) 

where 9 is the angle between polarization of light and 
the b axis. For an orthorhombic structure M x ^ M y . 
On the other hand, the main contribution to M which 
comes from atomic wave functions is independent of 9, 
so we can expand M in the vicinity of the Fermi energy 
and express matrix element in the form: 

M w , v (e)=M +y Xjy e (5) 

where j x y = dM x ^ y /de and the derivative is taken at the 
Fermi energy. Substituting this expression into Eq.(4) 
we obtain a qualitative description of the angular depen- 
dence of the probe signal G(9) and P{9): 

AR oc Aa cx M Q (~f x sin 2 (0) + ly cos 2 (9)) An (6) 

Here An is the number of photoexcited quasiparticles, as 
described in the Ref.[12]. The values of j x y can be dop- 
ing dependent and can easily have different sign, depend- 
ing strongly on w and the material's band structure in the 
vicinity of the resonance t — Ep ~ hu>±A. However, since 
the probe polarization anisotropy is a consequence of the 
anisotropy of the probe transition matrix elements, Eq. 

unfortunately does not give any direct information re- 
garding the anisotropy of the low-energy electronic gap 
structure. In Fig. 3 c) we plot the polarization depen- 
dence of photoinduced reflectivity amplitude given by Eq. 
@ using r y x /^/ y — 35/(— 10), which is shown to describe 
the main features of the data quite well, despite the fact 
that Eq.(^) is based on the oversimplified expansion of 
the matrix elements (Eq. ||) - neglecting higher order 
terms. 

An alternative model discussing carrier relaxation 
dynamics in HTSC has been recently proposed [^6| , 
whereby the main contribution to AR/R in the region 
at u> ~ 1.5 eV comes from spectral weight transfer in the 
real part of the optical conductivity a (to) from lo = to 
uj « A . Since Act is always positive this leads to negative 
changes in the real part of the dielectric function, inde- 
pendent of polarization. The change of sign in AR/R 
as a function of probe polarization which we observe in 



Y124 is clearly inconsistent with the model proposed in 
ref. @. 

In conclusion, the femtosecond relaxation dynamics in 
Y124 is found to be very similar to that reported previ- 
ously in Y123 [0. The results are also consistent with 
time-resolved terahertz spectroscopy measurements [^H) , 
which directly probed the recovery of the condensate af- 
ter photoexcitation, both agreeing well with the model 
Q. The observed probe polarization dependence and 
sign change of the transient signal below T c is described 
well by a model which considers the anisotropy of the 
probe transition matrix elements in Y124, and is not di- 
rectly related to the symmetry of the order parameter. 

The authors would also like to acknowledge A.Mironov 
of the Chemical Department Moscow State University for 
single crystal x-ray investigation of the samples. 



I. FIGURE CAPTIONS 

Figure 1: The photoinduced reflection AR/R from 
YiBa2Cu408 at different temperatures above and below 
T c as a function of time, measured a) with the polariza- 
tion of probe pulse in the direction parallel to the crystal 
axis a and b) parallel to (the direction of chains) the 
crystal axis b. The insert shows the data at 28K and 
74K already presented in a) on the logarithmic scale, so 
that two component decay can be easily observed. In 
c) we show weak oscillations of AR/R with frequency 
of 3.0(±0.3)THz observed above T c and attributed to a 
coherent phonon mode. 

Figure 2: The relaxation times tq as a function of 
temperature a) for YiBa2Cu40s, with the direction of 
polarization of probe pulse along a axis (open triangles) 
and b axis (squares), and b) for YBai.gSro^C^Os with 
the direction of polarization of probe pulse along b axis 
(squares). 

Figure 3: a) The temperature dependance of pho- 
toinduced amplitudes G(T) for YiBa2Cu40s, measured 
with the polarization of probe pulse along a axis (circles) 
and b axis (open triangles), and for YBai.6Sro.4Cu40s 
with the polarization along b axis (squares). The values 
of A(0) obtained from the fits are also shown, b) The 
photoinduced amplitudes P(T) for YiBa2Cu4C>8 and 
YBai.6Sro.4Cu40g as a function of temperature, mea- 
sured at the same polarizations as in a). We also report 
the obtained values of A p . c) The polarization depen- 
dance of photoinduced amplitudes G(T= 45 K, 9) and 
P(T= 45 K,9) for YiBa 2 Cu 4 8 , together with the po- 
larization dependence given by Eq. ^| (solid line). 

Figure 4: The four possible transitions, which can give 
rise to a photoinduced probe signal arising form a change 
in quasiparticle population An. If the transition proba- 
bilities do not cancel a non-zero photoinduced reflectivity 
transient is observed. The polarization anisotropy of the 
probe signal arises due to the transitions having different 
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probabilities in different directions relative to the crystal axes. Note: the scheme is drawn for T = OK. 
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